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Reverberation  Fluctuations  from  a  Smooth  Seafloor 


Steve  Stanic  and  Edgar  G.  Kennedy 


Abstract — High-frequency  shallow-water  reverberation  statis¬ 
tics  were  measured  from  a  smooth,  sandy,  featureless  seafloor. 
The  reverberation  statistics  are  presented  as  a  function  of  source 
frequency  (20-180  kHz),  grazing  angle  (30%  20°,  9.5C)  and  source 
beamwidths  (1.2=— 2.75').  Generally,  the  reverberation  statistics 
did  not  follow  a  Rayleigh  fading  model.  The  model  dependence 
of  the  reverberation  statistics  exhibited  a  complex  behavior  that 
ranged  from  near  Gaussian  to  beyond  log-normal.  The  results 
show  that  small  changes  in  the  source  frequency,  grazing  angles, 
and  beamwidths  caused  large  variations  in  the  model  dependence 
of  the  reverberation  statistics. 


I.  Introduction 

IGH-RESOLUTION  rough-surface  reverberation  mea¬ 
surements  have  shown  that  in  general,  Rayleigh  fading 
models  do  not  provide  an  adequate  description  of  rough- 
surface  clutter  statistics.  Rough  surface  reverberation  mea¬ 
surements  made  using  high-resolution  radars  and  sonars  have 
measured  scattered  envelope  amplitude  distributions  that  are 
described  by  Gaussian.  Rayleigh,  log-normal,  NVeibull  or  other 
more  complex  probability  density  functions  (PDF)  [l]-[9]. 
High-resolution  sonars  have  also  been  used  to  measure  seafloor 
reverberation  [10],  [11],  et  al.  Boehme  [12]  and  Stanic  et  al. 
[13]  have  shown  that  sea  floor  reverberation  statistics  also 
show  significant  departures  from  Rayleigh  fading  models.  In 
all  cases,  the  scattering  statistics  have  been  shown  to  be  highly 
dependent  on  sonar  frequency  and  beamwidth.  McDaniel 
has  related  Boehme's  and  Chotiros’  [II],  [12]  results  at  30 
kHz  to  the  measured  bottom  micro-roughness  spectrum  [14], 
Stanton  [15],  [16]  has  shown  that  the  shape  of  the  PDF 
can  provide  a  measure  of  the  relative  levels  of  coherent 
and  incoherent  scattering.  Stanton  has  also  related  the  PDF 
shapes  to  estimates  of  bottom  roughness,  correlation  lengths, 
and  bottom  microstructure.  These  results  have  .all  assumed 
that  the  seafloor  returns  were  not  contaminated  by  scattering 
from  subsurface  volume  scatterers.  Alexandrou  et  al.  [!7] 
have  compared  reverberation  PDF's  generated  by  point  scatter 
model  simulations  with  sea  beam  reverberation  data.  This 
comparison  lead  to  identifying  acoustic  signatures  to  be  used 
in  seafloor  classification. 

This  paper  presents  high-resolution  seafloor  reverberation 
statistics  as  a  function  of  frequency  (20—  1 80kHz),  grazing 
angle  (30°,  20°,  9.5°)  or  equivalent  range  intervals  (lengths 
of  insonified  bottom  area  at  the  — 3-dB  points),  and  system 
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and  by  the  Office  of  Naval  Research,  Program  Element  61 153N.  through  the 
NOARL  Defense  Research  Sciences  Program.  This  paper  is  NRL  Contribution 
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Fig.  I.  Schematic  of  insonified  seafloor  with  .V,  and  surface  and  volume 

scatterers. 

beamwidths  (1.2°-2.75°).  The  results  are  compared  to  various 
theoretical  HE^F  models  and  to  PDF’s  measured  for  a  coarse 
shell  bottom  during  the  Jacksonville  experiments  [13]. 

Fig.  1  shows  acoustic  energy  incident  on  a  surface  area 
A  that  contains  N,  surface  scatters.  In  addition,  the  acoustic 
energy  that  penetrates  into  the  sea  floor  is  scattered  by  a 
variety  of  individual  volume  scatters  A%.  For  the  monostatic 
backscattering  case,  the  signal  at  the  output  of  a  receiving 
sonar  sensor  is  given  by 

Vi  =  BiVi0  exp  [i(ur  -  2 kr{  +  0,-)] 

where  Vi0  is  the  backscattering  amplitude,  0;  is  the  scattering 
phase  of  the  signal  scattered  from  the  all  scatterer,  r;  is  the 
range  between  the  sonar  receiving  sensors  and  the  zth  scatterer, 
k  =  2-/A  is  the  acoustic  wave  number,  and  A  is  the  acoustic 
wave  length,  us  is  the  angular  frequency,  r  is  time,  and  5; 
is  a  system  constant  that  accounts  for  system  calibrations, 
propagation  losses,  and  other  system  related  factors.  If  the 
scatterers  are  independent,  then  the  total  reverberation  field 
due  to  the  Ns  and  Nv  scatterers  can  be  expressed  as 

jV 

F  =  J]5,oexp^t 
;=i 

where  $,•  =  usr  -  2 At,-  +  and  N  is  the  total  number  of 
scatterers. 

This  total  reverberation  field  will  in  general  fluctuate  around 
some  mean  value.  It  is  the  PDF  of  these  fluctuations  that  iden¬ 
tify  ihe  noise  models  against  which  target  detection  algorithms 
must  operate.  For  example,  a  processor  optimized  for  Rayleigh 
noise  would  perform  poorly  when  operating  against  a  non- 
Rayleigh  noise.  The  probability  of  false  alarm  (PFA),  defined 
as  one  minus  the  cumulative  distribution  function,  can  also  be 
used  to  characterize  reverberation  fluctuations.  The  PFA  is  a 
more  sensitive  measure  of  the  differences  in  the  tails  of  various 
PDF’s  than  the  cumulative  distribution  function. 

The  PFA  for  a  Rayleigh  PDF  is  given  by  [18] 

PFA(y)  =  cxp(-10y/1°). 
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Fis.  2.  Bo:tom  photograph  showing  ihe  smooth  featureless  seafloor. 


For  a  log  normal  distribution  the  PFA  is  given  by  (19] 


PFA(y)  =  /  exp 


<J\/2t7 


where  0  is  calculated  from 

wiocxP  [-((?/  ~  3)/a):]dy 


o  V  li 


a  is  defined  as  the  standard  deviation. 

The  PFA  for  a  Gaussian  PDF  can  be  expressed  as  [20] 


PFA(i/)  =  erfc 


y-m 

a 


where  m  is  the  signal  mean,  a  the  standard  deviation,  and 
erfc  the  error  function. 


II.  Experimental  Measurements 

A  series  of  high-resolution  bottom  reverberation  measure¬ 
ments  were  made  in  an  area  19  miles  south  of  Panama 
City,  FL.  The  seafloor  in  this  area  was  smooth,  with  only 
small  shell  fragments  spatially  distributed  in  the  sand.  Fig.  2 
is  a  photograph  of  the  sea  floor.  A  detailed  experimental 
description  of  the  measurement  area  is  siven  in  Stanic  et  al. 
[21]. 

Bottom  reverberation  measurements  were  made  using 
NRL’s  stable  acoustic  measurement  system  [22],  The  acoustic 
measurements  were  made  using  a  pair  of  nonlinear  parametric 
sources  and  a  two-dimensional  12-hydrophone  spatial  array. 
Source  beamwidth  as  a  function  of  frequency  for  the  250- 
kHz  and  450-kHz  sources  operating  parametrically  are  shown 
in  Table  I.  Data  from  the  receiving  hydrophone  located 
at  the  center  of  the  receiving  array  were  processed  using 
a  high-speed  computer  automated  measurement  and  control 
system  (CAM AC).  The  data  were  complex  base  banded  to 
5  kHz,  digitized  at  20  kHz  and  recorded  on  a  high-speed 
parallel  transfer  disk.  Individual  reverberation  envelopes  were 
calculated  and  the  results  recorded  on  optical  disks.  At  each 
grazing  angle,  the  source  frequency  was  varied  between  20 
and  180  kHz,  thus,  the  areas  insonified  by  the  450-kHz  source 
always  contained  the  areas  insonified  by  the  250-kHz  source. 


TABLE  II 

Horizontal  Ranch  Intervals  as  a  Function  of  Grazing 
Angle  and  Widest  Source  Beamwidth  usokHz  Source) 


Grazing 

angle 

Frequence 

kHz 

-5  dB 
beamwidth 

Horizontal 

range 

interval 

number 

Horizontal 
range 
interval 
length  (m) 

50  5 

90 

i — 

2.75= 

I 

12.4-13.9 

20° 

*90 

2.75° 

II 

19.4-22.5 

9.5 s 

90 

2.75s 

III 

39.5-53.2 

Table  II  shows  the  horizontal  range  intervals  for  the  widest 
source  beamwidth  (450  kHz7  at  each  grazing  angle. 

For  each  experimental  configuration,  reverberation  data 
were  taken  for  about  3  minutes  using  5-ms  long  continuous 
wave  (cw)  pulses  transmitted  at  1-s  intervals.  Since  the  re¬ 
verberation  data  were  taken  using  sonars  mounted  on  a  stable 
platform,  any  ping-to-ping  envelope  amplitude  variations  were 
caused  by  fluctuations  in  the  structure  of  the  water  column. 
These  fluctuations  caused  changes  in  the  size  and  position  of 
the  insonified  areas.  This  introduced  a  random  component  into 
the  reverberation  envelopes  that  increased  with  range. 

III.  Amplitude  Statistics 

The  data  taken  at  each  experiment  configuration  were  tested 
for  stationarity  using  the  Mann- Whitney  test  (23].  This  test  is  a 
nonparametric  test  that  can  be  used  to  determine  if  data  points 
form  part  of  a  locally  stationary  sequence  on  a  ping-to-ping 
basis.  For  each  experimental  configuration  3000  samples  were 
used  to  calculate  the  PDF's  of  the  reverberation  amplitude 
envelopes.  For  a  two  tailed  test  at  o  =  0.05,  i.e.,  95% 
confidence  level,  the  test  statistic  (Z„)  must  be  less  than  1.96 
in  order  for  the  data  points  to  be  considered  as  coming  from  the 
same  population  and  density  function.  Table  III  shows  the  first 
20  test  statistics  at  selected  frequencies  and  range  intervals. 
At  these  and  all  other  frequencies  and  range  intervals,  Zn  was 
less  than  1.96. 

The  reverberation  distribution  functions  were  normalized  so 
as  to  have  a  mean  of  one  and  plotted  as  the  probability  of 
false  alarms  (PFA).  These  results  are  compared  to  theoretical 
Rayleigh  and  log-normal  distributions  with  standard  deviations 
of  5.57  dB  [24J  and  to  a  Gaussian  distribution  with  a  mean 
of  zero  and  a  standard  deviation  of  one.  Fig.  3  shows  PFA 
curves  as  a  function  of  frequency  (180,  150,  110,  and  90 
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TABLE  III 

MannAVhitney  Test  Statistics  iz„)  for 
Selected  fr  equencies  and  Range  Intervals 


Frequency  ISO  kHz  Range 
Interval  I 

Frequency  90  kHzAVB  Range 
Interval  11 

Z\  =0.19647 

Z2  =  0.09996 

Z-2  =  0.10571 

Zs  =  0.0379.' 

Zj  =  0.11746 

Z*  =  0.09307 

Z5  =  0.34124 

Zi  =  0.01175 

Z.3  =  0.03170 

ZG  =  0.29365 

Zr  =  0.43544 

Z:  =  0.10571 

Z5  =  0.25552 

Zs  =  0.03524 

Z9  =  0.01034 

Z9  =  0.0469S 

Z,o  =  0.32036 

Zio  =  0.05222 

Zn  =  0.53427 

Zn  =  0.11746 

Zw  =  0.22405 

Zu  =  0.00397 

Zjj  =0.13511 

Z 1 3  =  0.03524 

ZM  =  0.51339 

Zu  =  0.14095 

Zi3  =  0.49291 

Zu  =  0.0469S 

Zie  =  0.27920 

Z  i  g  =0.17619 

Zit  =  0.52351 

Z  i  t  =  0.04695 

Zis  =  0.74109 

ZIS  =  0.02349 

Zis  =  0.72730 

Z19  =  0.2231S 

Frequency  60  kHz  Ranee  Interval 

I 

Frequency  20  kHz  Range  Interval 

111 

Zi  =  0.17274 

Z  i  =0.50212 

Z2  =  0.60459 

Z2  =  0.4S013 

Z3  =  0.19565 

Z3  =  0.11362 

Zi  =  0.05910 

Zi  =  0.35307 

Zi  =  0.35003 

Zj  =  0.26022 

Z6  =  0.12092 

Zo  =0.11362 

Zr  =  0.30229 

Zr  =  0.20591 

Z5  =  0.-30095 

Z*  =  0.00611 

Z9  =  0.33654 

Z9  =  0.125S3 

Zio  =  0.10364 

Zio  =  0.47524 

ZM  =  0.01727 

Zu  =  0.22113 

Zu  =  0.00564 

Zi'2  =  0.04276 

Z 1 3  =  0.06910 

Z13  =  0.0S1S5 

ZM  =  0.16410 

Zu  =  0.17470 

Z)5  =  0.00000 

Zn  =  0.44S36 

Z  is  =  0.25047 

Z  i  g  =  0.40194 

Z  i  t  =  0.04315 

Zir  =  0.23334 

Z 1 5  =  0.06046 

Zis  =  0.12339 

Z,s  =  0.11225 

Zj  o  =  0.20402 

kHz/\VB)  and  range  interval  for  the  data  taken  using  the 
450-kHz  source.  The  data  at  ISO,  150,  1 10,  and  90  kHz/WB 
and  at  all  range  intervals  showed  significant  departures  from 
Gaussian,  Rayleigh,  and  log-normal  models.  At  range  interval 
HI  the  data  at  150,  110,  and  90  kHzAVB  tended  to  follow 
a  PFA  other  than  log-normal.  At  range  interval  II  the  180 
and  150  kHz  data  tended  more  toward  Gaussian  than  did  the 
same  data  at  range  interval  I.  Only  at  90  kHzAVB  did  the  data 
show  a  sequential  model  dependence  with  increasing  range 
interval.The  PFA  curves  at  90  kHz/NB,  60,  40,  and  20  kHz 
arc  shown  in  Fig.  4.  These  data  were  taken  using  the  250-kHz 
parametric  source.  The  data  at  90  kHz/NB  and  range  interval 


Fig.  3.  Probability  of  false  alarms  as  a  function  of  frequency  and  range 
r~  intervals. 

r 


Fig.  4.  Probability  of  false  alarms  as  a  function  of  frequency  and  range 

imervals. 


II  tended  to  follow  a  Rayleigh  model.  At  range  intervals  1 
and  III  the  90-kHz/NB  data  tended  to  follow  PFAs  that  were 
characteristic  of  high-ground  clutter.  At  60,  40,  and  20  kHz  the 
data  at  all  range  intervals  showed  significant  departure  from 
ail  chosen  theoretical  PFA  models.  Only  at  60  kHz  and  for 
probability  of  false  alarms  less  than  -1  did  the  data  tend  to 
follow  a  Gaussian  curve.  These  results  are  also  summarized 
in  Table  IV. 


IV.  Beamwidth  Dependence 

Fig.  5  shows  the  PFA  dependence  on  system  beamwidth 
and  range  intervals  at  90  kHz.  The  data  were  taken  using 
the  450-kHz  and  250-kHz  source  transmitting  parametrically 
at  90  kHz  with  a  beamwidth  of  2.75°  (WB)  and  1.2°  (NB) 
respectively.  At  ranee  interval  I,  the  steepest  grazing  angle, 


TABLE  IV 

PFA  Mods-  Tendency  as  a  Function  of  Range  Interval  and  Frequency 


Range  Interval  I 


Model 

Gaussian 

Between 
Gauss¬ 
ian  and 
Rayleigh 

Rayleigh 

Between 

Rayleigh 

and 

log¬ 

normal 

log¬ 

normal 

Greater 

than 

log¬ 

normal 

Frcq- 

kHz 

ISO 

X 

150 

mm 

no 

H 

90WB 

X 

90NB 

X 

60 

X 

40 

X 

20 

X 

Range  Interval  II 


Model 

■ 

Gaussian 

Between 
Gauss¬ 
ian  and 
Rayleigh 

Rayleigh! 

Between 

Rayleigh 

and 

log¬ 

normal 

log¬ 

normal 

Greater 

lhan 

log¬ 

normal 

Frcq- 

kHz 

180 

X 

150 

X 

110 

X 

90WB 

X 

90NB 

X 

60 

X 

40 

x 

20 

X 

5 

u. 

S 

c 


6. 

8  '* 


a 

* 

& 

t 


a 

I- 

e 

I  0 
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Fis.  5.  Probabilitv  of  false  alarms  as  a  function  of  system  beamwidth  of 

90  kHz. 


Range  Interval  III 


Model 

Gaussian 

Between 
Gauss¬ 
ian  and 
Rayleigh 

Rayleigh 

Between 

Rayleigh 

and 

log¬ 

normal 

log¬ 

normal 

Greater 

than 

log¬ 

normal 

Frcq- 

kHz 

ISO 

X 

150 

X 

110 

X 

90  WB 

X 

90NB 

X 

60 

X 

40 

X 

20 

X 

the  90  kHz  (WB)  data  tend  to  lie  between  a  Gaussian  and 
log-normal  curve,  depending  on  the  value  of  PFA.  The  NB 
data  tended  to  follow  a  model  other  than  a  high-clutter  log¬ 
normal  model.  At  range  interval  II,  both  the  WB  and  NB  data 
tended  to  converge  on  the  Rayleigh  curve.  At  range  interval 
Hr,  both  the  WB  and  NB  data  showed  significant  departures 
from  a  log-normal  curve.  These  results  are  also  summarized 
in  Table  V. 


TABLE  V 

PFA  Model  Tendency  as  a  Function  of  Ranges 
Interval  and  Source  Beamwidth  at  w  ih t 


Model 

Gaussian 

Between 
Gauss¬ 
ian  and 
Rayleigh 

Rayleigh 

Between 

Reyleigh 

and 

log¬ 

normal 

log¬ 

normal 

Greater 

than 

log¬ 

normal 

Range 

Int. 

i 

ii 

IV 

WB 

NB 

WB 

NB 

NB. 

WB 

V.  Discussion 

The  results  presented  in  this  paper  show  that  reverberation 
fluctuations  measured  from  a  sandy  sea  floor  exhibit  large 
PFA  variations.  Even  so,  this  area  was  in  general  considered 
smooth  (non-resolvable  topographical  features).  Changes  in 
the  size  and  location  of  the  insonified  area  resulted  in  large 
variations  in  the  model  dependence  of  the  bottom  reverberation 
statistics  as  a  function  of  frequency,  range  interval,  and  system 
beamwidth. 

The  reverberation  statistics  measured  from  this  sandy  area 
had  different  model,  frequency,  and  range  interval  dependen- 
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cies.  than  the  reverberation  statistics  measured  from  the  shell 
covered  area  off  Jacksonville,  FL  [13].  At  the  Jacksonville 
site.^data  taken  at  ISO,  150,  and  110  kHz  and  at  range 
interval  l  tended  to  lie  between  Gaussian  and  Rayleigh  PFAs. 
The  Panama  City  data  taken  at  the  same  frequencies  and 
similar  range  intervals  had  a  different  model  dependence.  At 
frequencies  at  ISO,  150,  and  110  kHz,  the  model  dependence 
at  range  interval  II  was  between  Gaussian  and  Rayleigh. 

A  comparison  of  the  Panama  City  and  Jacksonville  PFA's 
at  90,  60,  40,  and  20  kHz  showed  that  the  Panama  City 
data  was  spread  between  Gaussian  and  log-normal,  while  the 
Jacksonville  data  was  in  general  spread  between  Rayleigh  and 
PFA's  that  were  characteristics  of  a  clutter  environment  that 
was  higher  than  that  of  a  log  normal  PFA.  For  the  Panama 
City  data,  only  one  PFA  curve  tended  to  follow  a  Rayleigh 
model:  90  kHz/XB  at  range  interval  II. 

For  the  Panama  City,  Jacksonville,  and  all  other  reported 
bottom  reverberation  data,  it  is  assumed  that  the  reverberation 
was  caused  by  scauerers  located  on  the  sea  floor  surface. 
Recently,  Boyle  and  Chotiros  [25]  have  shown  that  even 
at  angles  below  critical  a  considerable  amount  of  acoustic 
energy  can  penetrate  the  sea  floor  sediment.  The  critical  angle 
for  the  Jacksonville  data  was  26°.  Thus,  the  reverberation 
statistics  must  be  a  function  of  both  surface  and  volume 
scatterer  distributions.  These  scauerers  can  have  a  large  range 
of  sizes,  some  distributed  randomly  and  others  distributed  in 
all  types  of  uneven  and  patchy  configurations.  In  addition,  their 
scattering  facets  can  be  oriented  at  all  angles  with  respect  to 
the  maximum  response  axes  of  the  sonar  system.  Because  of 
this  complexity,  no  current  theory  can  completely  describe 
the  effects  of  the  various  scatterer  distributions  on  the  model 
dependence  of  the  reverberation  statistics.  The  authors  will 
attempt  to  address  some  of  these  questions  with  data  taken 
during  a  series  of  scattering  experiments  scheduled  for  the 
Gulf  of  Mexico  in  August  1993. 
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